Spectral structure and decompositions of optical states, and their applications 
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We discuss the spectral structure and decomposition of multi-photon states. Ordinarily 'multi- 
photon states' and 'Fock states' are regarded as synonymous. However, when the spectral degrees 
of freedom are included this is not the case, and the class of 'multi-photon' states is much broader 
than the class of 'Fock' states. We discuss the criteria for a state to be considered a Fock state. We 
then address the decomposition of general multi-photon states into bases of orthogonal eigenmodes, 
building on existing multi-mode theory, and introduce an occupation number representation that 
provides an elegant description of such states. This representation allows us to work in bases im- 
posed by experimental constraints, simplifying calculations in many situations. Finally we apply this 
technique to several example situations, which are highly relevant for state of the art experiments. 
These include Hong-Ou-Mandel interference, spectral filtering, finite bandwidth photo-detection, 
homodyne detection and the conditional preparation of Schrodinger Kitten and Fock states. Our 
techniques allow for very simple descriptions of each of these examples. 

PACS numbers: 42.50.-p,42.50.Dv,42.50.Ar 
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signal photon, but the spatial and spectral correlations 
between the photon pairs have a large impact on the state 
purity and thus on performance of quantum interfer- 
ences in networks @- Different methods for engineering 
the spectral properties of spontaneous parametric down- 
conversionprocesses for conditional single-photon prepa- 
ration 0, Q and entanglement based applications have 
been studied [l(| El • For bi-photon states the Schmidt 
decomposition has been proven to be a successful tool 
to analyze the properties the spatial and spectral de- 
grees of freedom of quantum fields [!, 0, [H, EH, E3 ■ The 
multi-mode character of multi-photon states is often ig- 
nored or may be irrelevant. Though for the generation 
of pulsed multi-photon squeezed states, spectral proper- 
ties and inter-mode correlations play an important role 
in squeezing optimization and for the complete charac- 
terization of realistic sources [lllEU]. Recent experiments 
now employ for the first time squeezed states in combi- 
nation with conditional state preparation to realize non- 
Gaussi an q uantum states with negative Wigner functions 
nam, |19( . These are crucial for further progress in CV 
quantum communication (20| and computation [21] appli- 
cations, but the theoretical modelling so far is frequently 
restricted to a single-mode representation. 

The need for more complex multi-photon states of light 
makes the accurate description of such states in an el- 
egant and compact representation an important issue. 
Most notably, in the context of quantum networks and 
conditional state preparation it is important to realize 
that different physical processes can impose a particu- 
lar basis choice, because they operate on only a particu- 
lar well-defined mode. For example, filtering and detec- 
tion processes are characterized by a particular spectral- 
temporal response. In this paper we analyze the inter- 
nal spectral structure of multi-photon states and study 
*rohde@physics. uq.edu.au URL: http : //www. physics .uq. edu. au/peopifctfiiafiipU'eations on quantum interference effects. We dis- 



I. INTRODUCTION 

Quantum communication and optical quantum com- 
putation rely on the controlled preparation and manip- 
ulation of specific, well-defined photonic states. In re- 
cent decades two main frameworks for the implementa- 
tion of quantum networks have been established: for sys- 
tems based on single-photon states the quantum infor- 
mation is encoded in mode properties (e.g. polarization, 
frequency or optical phases) of the optical fields, while 
it is normally assumed that each channel contains ex- 
actly one photon. Contrariwise, multi-photon states are 
typically employed in the context of continuous variable 
(CV) coding where the quantum observables are conju- 
gate quadratures which are directly linked to the photon 
numbers of the fields. In this case the description of the 
field properties is mostly reduced to a single-mode phase 
space representation [l[ ignoring spatio-spectral degrees 
of freedom. 

For single-photon states the impact of the spatio- 
spectral structure has been the subject of previous inves- 
tigation 0, [H, 0, [1[ . Experimentally there exist different 
approaches to realize single-photon states. If the emis- 
sion of single quantum systems is utilized, e.g. by em- 
ploying single atoms, molecules or quantum dot systems 
(for review see e.g. [6|), the suppression of any higher 
photon number states for each individual creation event 
avoids a complex internal structure. As a common al- 
ternative conditional single-photon states can be gener- 
ated by spontaneous parametric downconversion, which 
produces distinguishable photon pairs. In this case the 
detection of a trigger photon heralds the existance of a 
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cuss the properties, symmetry and decomposition of 
the spectral distribution functions of generalized states 
with higher photon number components. Additionally, 
we present techniques for performing multi-mode calcu- 
lations with such states. Previous work has introduced 
a multi-mode theory for the specific case of homo- 
dyne detection, where an optical field is decomposed into 
components orthogonal to and overlapping with some de- 
sired mode. We expand upon this notion with a gen- 
eralized eigenmode decomposition for arbitrary states. 
Our theoretical model provides us with great freedom 
in choosing the basis of the decomposition, while re- 
maining completely general. The flexibility in choosing 
the eigenmode decomposition basis makes it very suit- 
able for describing many physical systems of experimen- 
tal interest. In quantum networks the physical processes 
impose a decomposition in terms of the following ba- 
sis: One mode is defined by the response, and all oth- 
ers are moved to a different Hilbert space and may be 
traced out. Thus, the response of the filtering/detection 
can significantly change the nature of the state. This 
can have significant implications on, for example, con- 
ditional state preparation techniques, or detected char- 
acteristics linked to photon-number statistics, e.g. phase 
space representations reconstructed from homodyne mea- 
surements. Our approach reduces the complexity of cal- 
culations which involve higher photon number states 
and/or multiple modes and highlights the influence of 
the internal structure of pulsed photonic states. In par- 
ticular for multi-photon statesour analysis indicates the 
limitations of single-mode description where the frequen- 
cies and wavevectors of the involved fields are not taken 
into account. 

We begin by introducing a general representation for 
optical states, which includes both the photon number 
and spectral degrees of freedom in Section [TTJ We re- 
view the structure and normalization of such states and 
demonstrate that these states are far more general than 
the usual notion of 'Fock' states. We then analyze the 
relationship between generalized multi-photon states and 
Fock states, and give criteria for photonic states to be for- 
mally regarded as Fock states. We describe the impact on 
quantum interference effects. In Section Hill we introduce 
a decomposition of states into discrete, orthogonal 'eigen- 
modes', as well as an 'eigenmode occupation number' rep- 
resentation. This provides an elegant representation for 
general states, and is a useful tool when performing calcu- 
lations. For example, if we wish to understand a specific 
interaction between two arbitrary states, we proceed to 
decompose them into matched discrete eigenmodes. The 
overall evolution is then given by independently apply- 
ing the relevant interaction to each eigenmode. We thus 
take advantage of the fact that photons within a given 
eigenmode are by definition indistinguishable and can be 
treated in the idealized sense. In Section HVl we demon- 
strate the benefits of our technique by applying it to sev- 
eral examples of experimental relevance. The described 
decomposition techniques allow for very straightforward 



analysis of these model situations. While we treat the 
detection in a rather idealized way, it becomes apparent 
that already the generation of multi-photon states with 
internal spatio-spectral correlations prohibit a straight 
forward single-mode discription in many cases. Narrow- 
band filtering can only be applied in particular experi- 
ments to improve the purity of the quantum states. How- 
ever, our techniques could also prove very useful in the 
theoretical analysis of other quantum optical systems, 
in particular those subj ect to non-ideal effects such as 
mode- mismatch [2^, l24l [25l [26j j , photon distinguishabil- 
ity and imperfect photo-detection [27j • We conclude in 
Section El 

II. GENERALIZED REPRESENTATION OF 
MULTI-PHOTON STATES 

A. Theory 

We begin by discussing the representation of optical 
states in the spectral domain. Note that while we focus 
on the spectral degrees of freedom, the discussed repre- 
sentations can easily be generalized to other degrees of 
freedom, such as the transverse spatial ones. 

Including the spectral degrees of freedom, a pure n- 
photon state can be expressed in the form 

\ip n ) oc J $(ui,...,w n )tf(u 1 )...a\w n )dw\0), (1) 

where aJ(u>) is the single frequency photonic creation op- 
erator at frequency uj. i\> is a continuous function which 
characterizes the spectral distribution of the state. We 
refer to this as the spectral distribution function (SDF). 
An arbitrary pure state can then be expressed in the form 

\f) cx^c'„ / V'™(wi,...,w„)a t (wi)...a t (o;„)dd;|0). 

n J 

(2) 

We now turn our attention to the normalization of such 
states. Our discussion closely follows the work of Ou [HI]. 
Let us assume the SDF is initially normalized according 
to 

J ■ ■ ■ , uj n )\ 2 duj = 1. (3) 

For a single photon state this is equivalent to requiring 
that the state itself be normalized, i.e. (ip\ip) = 1. How- 
ever, for states with higher photon number this is not nec- 
essarily the case (i.e. in general (ip\ ip) ^ 1), as we now 
discuss. For the n-photon state \ip n ), we introduce the 
normalization term, Af n , defined such that (tp n \ipn) = 1. 
Thus, a normalized multi-photon state is of the form 

\tp) = —= / ip n (u)i, . . . Lu n )a\wi) . ..a){u n ) dw|0), 

(4) 
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where ^2 n |c„| 2 = 1. Following Ref. (28|, the normaliza- 
tion terms are given by 

Afn = \ ^n{^l, ■ ■ ■ , W n )* ^2 1pn(P[u>l, ■ ■ ■ ,W n ])dw, 

Pes„ 

(5) 

where P represents a permutation of the indices of ?/>, 
and we sum over all possible permutations. This expres- 
sion arises directly by expanding {ip n \ ipn}- This expan- 
sion consists of n! terms, containing all possible combi- 
nations of n-fold products of delta functions of the form 
(uji\u)j} = 8{uji — u!j) V i, j = 1 . . . n. This gives rise to the 
summation over permutations of the indices. Since there 
are n\ terms in the summation, the factor Af n ranges be- 
tween 1 and n\, and directly gives us a measure of the 
permutation symmetry of if). 

To provide some intuition into the behavior of Af, let 
us consider two simple examples. First, suppose the SDF 
can be expressed in the form 

^(wi,...,w„) = <t>{ui)...(j){u n ). (6) 

In this case ip exhibits full symmetry under per- 
mutation of the indices, i.e. ip(P[u>i, . . . ,u n }) = 
i/>(u>i, . . . , u n ) V P G S n , and it follows from Eq. [5] that 
Af = nl. It is shown in Appendix [Al that Af = nl if and 
only if tp exhibits full permutation symmetry. Note that 
full permutation symmetry does not uniquely correspond 
to the class of states of the form shown in Eq. [6l i.e. in- 
distinguishable separable states. For example, if we take 
the SDF to be a linear combination of two distinct such 
functions, e.g. ^(wi,cji) = a(f>(u>i)(j)(oj2) + P^p(oJi)ip(oj2), 
it will also exhibit full permutation symmetry, but is non- 
separable. 

Next, suppose the multi-photon state consists of 
temporally-separated photons, sufficiently separated 
such that every photon is distinguishable from the others. 
This essentially corresponds to n photons in n distinct 
modes, as opposed to n photons in one mode as per the 
previous example. Now the SDF can be expressed in the 
form 

w n ) = e^(«i)e-^ T ...e-^ T ^) 

n n 

= n e ~ fe " fcT ^)=n^M' ( 7 ) 

k=l k=l 

where r is the temporal separation between photons and 
is larger than the temporal bandwidth of the individual 
photons, i.e. r ^> 1/Aw. Thus, 

J ipi(u)*ipj(u) « 5 it j. (8) 

In this case it is obvious that the state exhibits no per- 
mutation symmetry, since the state decomposes such that 
each photon effectively has zero overlap with all others. 
Thus, the only surviving term in the expansion is for 
P = 1, giving Af = 1. 



Note that in the previous two examples the Af nor- 
malization factor is consistent with what we expect from 
standard theory. In the former case we have n excitations 
of a single mode, in which case standard theory tells us 
that \tp) — -^=(a^)"|0). On the other hand, in the latter 
case we effectively have n distinct modes, each with a 
single excitation, for which \ip) = a\a 2 ■ ■ ■ ^n|0). 

This notion generalizes very simply to the multi-mode 
case. Consider the m-partite system with rii photons in 
the ? th mode. Such a state can be expressed in the form 

xa+(u,«) . . . fit«>) . . . aUuM) . . . at>M) dtf|0). 

(9) 

where (df' denotes the frequency of the i th photon in 
mode j, and aj(w) denotes the single frequency creation 
operator of mode i. In this case it follows that the nor- 
malization factor Af is given by 




where Pi is a permutation over the just the indices of 
mode i, and again we sum over all permutations. It can 
easily be seen that in the multimode case Af varies be- 
tween 1 and ni\ri2\ ■ ■ ■ n m \. 



B. What is a Fock state? 

It is evident that there is an entire class of states con- 
sisting of n photons. This raises the question as to what 
subset of these states is equivalent to the usual notion of 
Fock states. We answer this question by examining the 
algebraic structure of Fock states. Formally, the n photon 
Fock state is defined as 

|n) = ^(«T|0). (11) 

To maintain a consistent algebraic structure we require 
that states defined in the more generalized spectral rep- 
resentation exhibit an analogous structure, 




where is a generalized mode creation operator, which 
creates a photon characterized by SDF 4>{u>), 

1+ = / ' $(u)tf{w)du>. (13) 

It is clear that the class of states corresponding to the 

Fock states is defined by the constraint on the multi- 
photon SDF, 

ip(uJi, ...,w n ) = cj>{wij(j}{u)2) ■ ■ ■ (p(u„), (14) 
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where 4>{oj) is arbitrary. Thus, when the SDF of an n- 
photon state is factorizable as an n-fold product of iden- 
tical functions, it is structurally equivalent to a Fock 
state. Henceforth we will refer to states of this form 
as Fock states. Note that for Fock states N = nl, giv- 
ing rise to the required normalization factor of \j\fn\ in 
Eq. |T2] Importantly, although Fock states necessarily sat- 
isfy M = nl, states satisfying Af — nl are not necessarily 
Fock states. An example of this are entangled states of 
the form ip(uii,u!i) — a(f>(oji)(/)(uj2) + /3 < p( w i)v( w 2) ) that 
were discussed earlier. Such states are fully permutation 
symmetric but are not Fock states. 

Often spectral distinguishability of multi-photon states 
is ignored, and an n-photon Fock state and a state with 
n photons are regarded as being equivalent. As we have 
discussed, Fock states only correspond to a subset of the 
later more general class of states. It is not immediately 
obvious what implications this has in an experimental 
situation. An example of this question was studied in 
Ref. [28[ , and temporal distinguishability in multi- photon 
states has been experimentally demonstrated in Ref. [2^ | . 

To illustrate this, and provide some intuition into the 
importance of the distinction between Fock states and 
more general multi-photon states, we now discuss a very 
simple example situation where the behavior of an ex- 
periment differs significantly, depending on the spectral 
composition of multi-photon states. 

Consider an experiment where we interfere two iden- 
tical copies of a two-photon state on a 50/50 beam- 
splitter. We will consider two opposing limits. First, we 
employ Fock states, as defined in the previous discus- 
sion. In this case each of the incident two-photon states 
are of the spectral form i/j(uji, w 2 ) — 0(wi)0(w2)- Sec- 
ond, we employ two-photon states in which the con- 
stituent photons are well temporally separated. In this 
case, ip(u)i,u>2) = e~ lUlT </>(wi)</>(w2), where r is much 
larger than the photons' temporal bandwidth. 

Let us consider the probability of detecting four pho- 
tons at output mode A, P^a- In the former case the 
derivation proceeds in the usual manner. 



IV'in) 

IVw-t) 



= U\ll>in) 

(fit) 4 (fit) 2 ^) 2 (St) 4 



|o) 



||4,0> - ||2,2) + a/|o,4), 



(15) 



where U is the 50/50 beamsplitter operation and we 
assume the phase-asymmetric beamsplitter convention. 
Therefore P 4A = 3/8. 

In the later case, temporally disjoint photons do not 
interact with one another at the beamsplitter. We can 
therefore treat the evolution as two independent interfer- 
ences of single photons at the beamsplitter. For a single 



such event we have 

IV'in) 

IVw) 



= at&t|o) 

'(fit) 2 



(St)* 
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= >'°>->' 2 >- (16) 

Thus, for each event the probability of detecting two pho- 
tons at any one of the outputs is 1/2, as expected from the 
well known Hong-Ou-Mandel (HOM) effect [13]. Given 
two independent events of this type, the total probabil- 
ity of detecting four photons at output A is P4A = 1/4. 

Although in both cases we are interacting two indis- 
tinguishable states, it is evident that the nature of the 
quantum interference differs significantly depending on 
their internal structure. 

These two limiting cases are easily calculated from 
standard theory. However, we lack a description of the 
more general case where the two-photon states exhibit 
partial internal distinguishability. We now consider a 
more general case. Suppose our incident two-photon 
states are separable and of the form 



I 02 



A l A U°)- 



(17) 



where ipi (lo) and ip 2 {u>) are arbitrary and may exhibit any 
degree of distinguishability. It can easily be seen that the 
normalization factor A2 is given by 

A/2 = ( / ipi(uj)*Lp 1 (uj)duj} I / tp 2 (u)* ip 2 (uj) du j 



+ \J fii^)* ip 2 {uj) dco j I J ip 2 (uj)* ipi(uj) duj 
= 1 + 7, (18) 



where 
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tpi(uj)* ip 2 (uj) du> 



(19) 



characterizes the degree of photon distinguishability 
within the two-photon state. We have two identical copies 
of this state, so the incident state is of the form 



|0n 



1 



(20) 



where A and B denote distinct spatial modes. We then 
pass this state through a 50/50 beamsplitter and con- 
dition upon detecting all four photons at the A out- 
put. When applying the beamsplitter transformation we 
assume that corresponding spectral modes from spatial 
modes A and B, match one another at the beamsplitter 
(we make this assumption throughout this paper). It can 
now easily be seen that the conditional output state is 



IV'cond) 



1 



(21) 
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The probability of detecting this state is simply given by 

PiA = (VWidl^cond)- Thus, 



PiA = 



A/~4 



(22) 



where A/4 is the normalization factor of the four-photon 
state A^A^A^A^JO). Notice that the interference is 
completely characterized by the relationship between dif- 
ferent permutation symmetry factors. It can easily be 
shown (see Appendix [C] for proof) that in this case 
A/4 = 4(1 + 7+ 7 2 )- Thus we have 



Pa 



1 I 1 



(23) 



As expected, this is consistent with the limiting cases con- 
sidered previously. When we have incident Fock states, 
7 = 1 and P4A — 3/8, whereas for states in which the 
constituent photons are distinguishable we have 7 = 
and P 4A = 1/4. 



where X^J^I 2 = 1' an d the discrete set of functions 
{£i(w)} are orthonormal under the inner product, 



(25) 



In the present discussion the exact choice of eigenfunc- 
tions will typically be unimportant. However, we will 
extensively make use of one property: with an appro- 
priate change of basis transformation we can choose an 
arbitrary normalizable function to belong to the dis- 
crete set of basis functions. This derives from the Gram- 
Schmidt orthogonalization procedure and is proven in 
Appendix [B] 

It follows that a single photon state with arbitrary SDF 
can always be expressed as a superposition of single pho- 
ton states with orthogonal SDF's, 



i 

= $>ijj0>. (26) 



III. EIGENMODE DECOMPOSITION OF 
MULTI-PHOTON STATES 

A. Motivation 

From the previous section it is evident that the internal 
spectral structure of multi-photon states plays an impor- 
tant role in quantum interference effects. The calculation 
of the previous example also illustrates that it becomes 
increasingly difficult to model generalized settings with 
higher photon numbers, since we must analytically ma- 
nipulate complicated multi-fold integrals with increasing 
numbers of terms. In this section we discuss an alterna- 
tive approach for representing multi-photon states in a 
discrete basis of eigenfunctions, which provides an ele- 
gant and compact description for such states and simpli- 
fies many calculations. We will demonstrate its benefits 
in the next section by analyzing several example situa- 
tions. Our method lends itself well to numerical analysis, 
as opposed to analytic manipulation of states in inte- 
gral form. This approach is an extension of the multi- 
mode description of optical states, which has been used 
in the evaluation of homodyne detection [22| , and is also 
closely related to the well-known Schmidt decomposition 
for two-mode states. 



B. Theory 

Any well-behaved complex function (in the present dis- 
cussion a SDF) can always be decomposed in terms of a 
discrete basis of orthonormal functions (a well known ex- 
ample is the basis of Hermite functions), 



(24) 



Because the basis functions {£i(w)} are orthonormal, 
(0| j4| |0) = Sij. Thus, the modes characterized by 

the different mode creation operators A|. are orthonor- 
mal. We refer to these as 'eigenmodes'. The coefficients 
{Xi} follow trivially from the overlap Aj = (0|j4^^4v|0), 



Xi 



£j(w)*V>(w) dui. 



(27) 



This notion generalizes to states with higher photon num- 
ber in a straightforward manner, 

\ip n ) = J V>(wi,...,w„)a t (wi)...a t (w„)dw|0) 
= K u ...,i n 

ii<.-.<i n 

x / . . . €i n (w n )a! (wi) . . .a t (w„)dw|0) 



(28) 



Notice that we insert the summation condition i\ < • • • < 
i n to remove potential ambiguity in the decomposition. 
^ii....,i„ is to be interpreted as the coefficient of the term 
with a photon in each of the modes ii, . . . , i n . Therefore 
A's with different permutations of the same indices refer 
to the same eigenmode. For example, in the two photon 
state the coefficients Ai^ and A2,i are both associated 
with a single photon in each of the modes £1 and £2- Thus, 
to avoid this double counting we assume the summation 
does not run over permutations of the same indices. Now 
the coefficients of the decomposition are uniquely defined 
as 



Xi 



, u) n ) duj. 



(29) 
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We can reformalize this representation slightly as fol- 
lows. Consider the single photon case. We define = 

Ac. |0) to be the single photon state in the eigenmode 
characterized by £i. Because the different eigenmodes are 
orthonormal, we can introduce an occupation number 
representation, 

|1>6 =|0> £l |0> & ...|l) 6 |0) £l+1 .... (30) 

Thus, an arbitrary single photon state can always be ex- 
pressed in the form 

|^i)=^A i |0) £l |0) ?2 ...|%|0) £i+1 .... (31) 

i 

Similarly, an arbitrary n photon state can be expressed 
in the eigenmode occupation number representation as 

|-0n) = E V n iW- ■ ■ ■Ai 1 ,..., i Jni)£ 1 |n 2 >£ 2 • • • 

ii<---<i n 

= E A U..,i>i>«>2>£ 2 ..., (32) 

i 1 <---<i n 

where rij is the number of indices . . . , i n } equal to j, 
i.e. the population of the respective eigenmode. This ef- 
fectively represents an arbitrary state in terms of a super- 
position of distinct Fock states. This change in notation 
may seem rather pointless, however this is exactly what 
we ordinarily do when we use regular Fock notation - we 
abstract away the specific form of the underlying wave- 
function to simplify treatment. This notational change 
turns out to be very useful, as we will see in the exam- 
ples in the next section. 

Next we discuss a slight variation of the eigenmode 
occupation number representation that simplifies certain 
types of calculation. This will be particularly useful in 
some of the examples we discuss shortly. In many calcu- 
lations we do not explicitly need to consider all of the 
modes in a given decomposition. Instead, there might be 
just one mode which is involved in a given interaction, 
and others are of no direct interest. In such a situation 
it would be particularly desirable to decompose a state 
into just two components: an arbitrary mode that is of 
interest (e.g. this mode might be involved in a particu- 
lar interaction), and collectively 'other' modes that are 
not directly involved in our calculations. For example, 
suppose we measure a state using a photo-detector that 
responds to a finite spectral range. By defining one of the 
eigenmodes to be the spectral response of the detector, 
we can simply trace out the remaining modes to under- 
stand the behavior of the detection process - the exact 
form of the rest of the decomposition is not important. 
This is one of the examples we will consider in detail 
shortly. 

First consider the single photon case. We represent the 
single photon SDF as, 

V>(w) = \i<f>(u) + (33) 



where <j)(u>) represents some 'desired' mode, and 4>{lo) col- 
lectively represents the component ofijj(u)) that does not 
overlap with <j>{u)). Now the mode creation operator can 
be recxpressed, 

4 = Ar4 + A it. (34) 

where 

Ai = J (f>(w)*il)(u)du), 

A = Vl- I^il 2 . (35) 

Moving to the occupation number representation this can 
be written 

|Vi> = Ai|l)*|O>0 + Ao|O> |l>0. (36) 

Of course this can also be easily generalized to states with 
higher photon number. For example, the decomposition 
of an arbitrary n photon state in terms of Fock states in 
mode (j> will take the form, 

n 

Thus far we have specifically focussed on the single- 
mode case. However, these ideas easily generalize to the 
multi-mode case. Specifically, a state distributed across 
modes A, B, C, . . . , with rii photons in mode i will take 
the form, 

|V»*,n fl ,...) = E X i A { A --- A l, 
i A .i B ,... " A 

xB\. B ...B\. B x...|0). (38) 

*1 l "B 

This representation for multi-photon states is closely re- 
lated to the Schmidt decomposition for biphoton states. 
Formally, a Schmidt decomposition corresponds to the 
tia = 1, Ub = 1 case where we choose the eigenbasis {&} 
so as to diagonalize the A matrix. Unlike the Schmidt de- 
composition, the representation described is completely 
general and holds for states with arbitrary photon num- 
ber. 



IV. EXPERIMENTAL EXAMPLES 

The eigenmode decomposition and occupation num- 
ber representation provide an elegant tool for perform- 
ing calculations with generalized multi-photon states. We 
now present several examples of experimental relevance, 
which illustrate the simplicity of modeling such systems 
using this representation. First we review the well known 
Hong-Ou-Mandcl effect, which provides a simple exam- 
ple scenario for our techniques. Second, we model spec- 
tral filtering and finite bandwidth photo-detection, which 
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follows very simply from an occupation number represen- 
tation, and is important in understanding the behavior of 
realistic photo-detectors. In the framework of this paper 
we simplify the calculations by assuming that the detec- 
tor projects onto a particular, well-defined mode |£o)- I n 
practise, realistic photo-detectors may respond to multi- 
ple modes but not provide information as to which one 
is detected. For example, a detector may respond to all 
times within some window, project onto a particular tem- 
poral mode within that window, but not tell us which 
one. In this instance we will clearly be left with a mix- 
ture of states of different possible detection outcomes. 
However, if we arc interested how internal correlations of 
multi-photon states impact the general performance an 
quantum optical system - derived in a single-mode pic- 
ture - an idealized detection provides us with an upper 
bound for the purity of the states. In the third section 
we discus homodyne detection in an effective single-mode 
description where we introduce the local oscillator mode 
as the relevant eigenmode. The fourth and fifth section 
finally treats two distinct conditional state preparation 
scenarios: the preparation of Schrodinger Kitten states 
via photon subtraction, and conditional Fock state prepa- 
ration via non-degenerate parametric down-conversion. 
These two examples illustrate that the achievalbe purity 
depends on the internal structure of the initial states as 
well as on the specific setups. 



A. Hong-Ou-Mandel interference 

We begin our discussion by considering the well-known 
HOM effect [30j ]. Although this effect is well understood 
and extremely simple to derive using various other ap- 
proaches (cf. Ref. [23(] ) , for illustration we start by red- 
eriving this result using the techniques discussed. After 
studying the case of independently produced photons, we 
generalize to the case of spectrally entangled photons. 

We begin with a two mode input state, with a single 
photon in each mode. We initially assume the state to 
be separable, such that the photons are independently 
characterized by SDFs 4>(uj) and f{to) respectively. The 
incident state can therefore be expressed in the form 



ID 



(39) 



where A and B denote distinct spatial modes. Shifting to 
the occupation number representation and choosing one 
of the eigenmodes to be 4>{oj) we obtain 



|V>)in = |1}«^4(Ai|1) ,b + A o |% b ) 



Ai|l) 



h.A 



I) 



Aoii; 



(40) 



Here we have used shorthand notation and not explic- 
itly written out all the vacuum terms. We now apply 
the 50/50 beamsplitter operation U . Note that because 
the different eigenmodes are orthogonal, and within each 
eigenmode photons are by definition indistinguishable, U 



acts on each eigenmode independently and in the ideal- 
ized sense. Thus we have the property, 



[/■(|mi,ni) ei |m2,n 2 )| 2 ...) 
= (U\m 1 ,n 1 )) il (U\m 1 ,n 2 ))^ 2 

Using this we obtain 

IV^out = U\lp)in 



(41) 



A 1 l7|l,l)* + Ao(l7|l,0))*(r7|0 J l>)$ 
^(|l,O) + |O,l)) (|l,O)-|O,l))^, (42) 



= (|2,0)-|0,2)k 



where we have used the shorthand \m,n)^ — 
\m)</>,A \n)$,B- From this expression we can directly read 
off the two-photon coincidence probability P c , 



lAol 2 



1 1 

2 ~ 2 



(43) 



This is the usual result - for completely indistinguishable 
photons the integral term approaches unity and P c = 0, 
whereas for completely distinguishable photons the inte- 
gral term approaches zero and P c = 1/2. 

This analysis effectively describes the situation where 
we interact two independently produced photons. Next 
we consider the more general case of spectrally entan- 
gled photons. This arises, for example, when using photon 
pairs produced through spontaneous parametric down- 
conversion. By expressing an arbitrary two mode state 
as 



(44) 



we can apply the 50/50 beamsplitter operation to obtain 
2 

•^A«[(i!) 2 -(£]) 2 ]|0>. (45) 



L>|^in) 

lY, X ^ A l+Bl)(A]~B])\0) 



Post-selecting upon coincidence events leaves 

l^cond) = ^^M]^ - A\B])\Q) 



Thus, the coincidence probability is 



(46) 



(47) 
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We can now immediately understand the behavior of 
HOM interference in the general case. First, suppose the 
input state is separable and the photons are indistin- 
guishable. Now it is necessarily possible to choose an 
eigenbasis such that An = 1, otherwise 0. In this case 
clearly P c — 0, as expected. Next consider the other ex- 
treme, where the two photons are completely distinguish- 
able. Suppose that photon A resides in eigenmode £i and 
photon B in eigenmode £2- Now A12 = 1, otherwise 0. 
Clearly now P c — 1/2, as expected. More generally we 
see that the requirement for perfect HOM interference is 
I Ay — Aj-i| = V This is equivalent to requiring that 
the incident two-mode state be exchange symmetric. In 
other words, the two-photon joint-SDF tp(cji,u>2) must 
be symmetric about the u)\ — L02 axis. Note that this is 
much broader than the class of indistinguishable separa- 
ble states, and includes highly spectrally entangled ones. 
This explains, for example, why HOM interference works 
with photon pairs produced through spontaneous down- 
conversion, despite the inherent spectral entanglement in 
that case (see also Ref. (3l|). 



B. Spectral filtering and finite bandwidth 
photo-detection 

We now consider the situation where we have an n- 
photon state to which we apply a spectral filter. Suppose 
the filter transmits frequencies within some finite range, 
and completely reflects any spectral components outside 
this range. Thus, the filter has an ideal rectangular spec- 
tral response function that we will label (j>(ui). For sim- 
plicity we will assume the incident state is an arbitrary 
Fock state. Note that this technique is by no means lim- 
ited to dealing with Fock states, it is just a simplifying 
assumption for the sake of example. The incident state 
can be expressed, 




where we have used an occupation number representation 
and chosen one of the eigenmodes to be (j>(u>). Next, by 
definition our filter reflects all modes orthogonal to (/>(uj). 
Thus, we trace them out to obtain the transmitted state, 

/3 = £("Wa 2 (— >|z)^U, (49) 
i=o 

which is a mixture of different photon number terms. 



Hence we find that in contrast to previous anaysis for 
bi-partite systems we find that in general spectral filter- 
ing does not improve the purity of the system but in- 
troduces mixing in the number degree of freedom. This 
is a generic property of filtering and imposes limitations 
on the applicability of filtering techniques. The notable 
exception to this is in all coincidence- type experiments, 
where photon number conservation conditions eliminate 
mixing in the number degree of freedom by effectively 
projecting us back onto a state with known photon num- 
ber. For example, present-day in-principle demonstra- 
tions of linear optics quantum computing (LOQC) [HJ 
operate in coincidence, enabling very narrowband filter- 
ing to post-select on highly overlapping spatio-temporal 
modes (see, for example, Refs. [33|, [3J, |35|). However, 
because these experiments operate in coincidence they 
are inherently destructive. Scalable implementations typ- 
ically require heralded quantum gates, in which case such 
filtering is ruled out due to the described mixing effect. In 
this context preserving purity of the states is extremely 
difficult. This poses one of the major challenges to the 
implementation of scalable LOQC. 

Next consider the case of finite bandwidth photo- 
detection. A finite bandwidth photo-detector can be 
modeled by considering an ideal detector preceded by 
an appropriate spectral filter. Thus, the state shown in 
Eq. 25] represents the state observed by the ideal detec- 
tor. The photon number statistics at the detector will 
simply be given by 

P m =(l\\x 2m \l {n - m \ (50) 

compared to the ideal case where P m = 5 m ,n- See 
Ref. [27| for a different treatment of finite bandwidth 
photo-detection. 

C. Homodyne detection 

Next we consider homodyne detection in an effective 
single-mode description which imposes mode-mismatch 
if we are not able or allowed to chose the local oscil- 
lator mode as underlying spatio-spectral mode. Exper- 
imentally we encounter this situation if we consider a 
quantum network with several - possibly correlated - 
signal modes which exhibit different spatio-spectral char- 
acteristics. In such systems we have to decide on a refer- 
ence mode for all channels in the network. The derivation 
of homodyne detection in presence of mode mismatch is 
well-understood (cf . Ref. [22j, |36[ ) , but we include this ex- 
ample to illustrate how our mode decomposition method 
relates to known results and to emphasize that spectral- 
spatial filtering can be hidden in the detection. 

Homodyne detection is widely used in quantum optics 
experiments and forms the basis of Optical Homodyne 
Tomography (OHT) [37]], which is used to tomographi- 
cally reconstruct the Wigner function of unknown states. 
The measurement proceeds by interacting an unknown 
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state with a coherent probe beam. Here an indirect fil- 
tering process takes place. Specifically, only the compo- 
nent of the incident state overlapping with the coherent 
probe mode will contribute to the homodyne statistics. 
As in the previous example, this effectively traces out 
the non-overlapping components. We now consider this 
in detail. 

For simplicity, let us consider a superposition of Fock 
states of the form 



(51) 



that we wish to measure. We will assume the coherent 
probe |q;) is formally a superposition of Fock states, char- 
acterized by mode function <j>(u>). We begin by reexpress- 
ing our state |^>) in the eigenmode occupation number 
representation, where we choose one of the eigenmodes 
to be (j>(u>) (the others are irrelevant for the calculation). 
Thus, 



= E C "E ^n,i\i)^\n - 



(52) 



n i—0 

and the corresponding density matrix is simply 

p = \i>){i>\ 

= Cn*c n >\ n ,i*Xn',i'*\i)4>\n - i)^(n' - i'l^ii'lc/,. 

(53) 

As discussed, the components of the incident state which 
do not overlap with the coherent probe do dot contribute 
to the homodyne statistics and are therefore effectively 
traced out. Thus, the observed density matrix is given by 

Pobs = tr^(/5) 

= Cn*C n > \n,i*^n',i'*5n-i,n>-i'\i){i'\- (54) 

Consider the behavior of this observed state in the limit- 
ing cases. First consider the ideal case, when the incident 
state is perfectly mode-matched to the coherent probe. 
In this case A 



>■:) 



Si j and it follows that 



Pobs 



= ^c„*c„»(n'| = 



(55) 



as expected. In the other extreme the incident state is 
completely mismatched from the coherent probe. Now 
X n ,i = &i.o an d the observed state is given by 



Pobs 



n,n f 



-|0)(0| = |0)(0|. 



(56) 



This is also intuitively expected - in the limit where 
the incident state is completely mismatched from the 
probe, the detection process just observes a vacuum state. 
The intermediate cases may be readily calculated given 



knowledge of the A's, which are easily derived from Eq.[29l 
for known ip. 

Notice that in this calculation we did not need to ex- 
plicitly model the homodyning process. Instead, it is ab- 
stracted away by decomposing the state into components 
which can be treated in the ideal sense, and which can be 
discarded. This is one of the main advantages of our ap- 
proach, and highlights the underlying physics inherent in 
the structure of multi-mode optical fields. If we want to 
study the properties of multi-mode states independent 
of their spatio-spectral character in an effective single- 
mode model we have to adopt this strategy to eliminate 
unobserved modes. 



D. Preparation of Schrodinger Kitten states 

We now turn our attention to a more specific exam- 
ple, a single-mode scheme for preparing Schrodinger Kit- 
ten states through photon subtraction from a pulsed 
single-mode squeezed state. This has been experimen- 
tally demonstrated in Ref. [13]. In this experiment a 
squeezed state is prepared via degenerate parametric 
down-conversion. A photon is non-deterministically sub- 
tracted from this state by inserting a low reflectivity 
beamsplitter and conditioning on detecting a single pho- 
ton in the reflected mode. When this post-selection pro- 
cess succeeds, we ought to have a photon-subtracted 
squeezed state in the transmitted mode, which corre- 
sponds to the desired Kitten state. We now examine the 
effect of the SDF on the preparation of such states. 

For a pulsed system we model the Hamiltonian of the 
parametric down-conversion process as 

H(t) — X / / -0(wi,w 2 )a t (wi)d t (w 2 )dwidw2 + H.c, 



where ip(w\,U32) is determined by the crystal properties, 
phase matching conditions, and pump pulse [111 ]. As we 
show in Appendix [D] the output state leaving the para- 
metric down-conversion crystal is given by 



E 

71 = 



A" 



ij){u>i, uj2)a) (wi)(r (W2) dtoi dix>2 



(58) 



Rewriting this in terms of an eigenmode decomposition 
we obtain 



E 

71 = 



A" 



EM14 



|0>. 



(59) 



Next we pass this state through a low-reflectivity beam- 
splitter, 



EME^A/i^I+y^t; 

71— ij 

A , A , \ -1 n 

y/T=jA]+y/TjB})\ |0), (60) 
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where r\ is the beamsplitter reflectivity. We then post- 
select upon detecting a single photon in the reflected 
mode B. In practise [171 ] . very narrowband filtering is em- 
ployed prior to photo-detection. Suppose the filter has an 
ideal top-hat response. Furthermore, we will assume that 
our eigenbasis decomposition is in terms of a basis where 
the first basis function £o is the filter response function. 
In the limit of small beamsplitter reflectivity, 77 <C 1, the 
conditioned transmitted state can be approximated as, 



E 

71=0 



A" 



(n-1)! 



EM 



E^ 



(61) 

In general this state is not a superposition of Fock states, 
and therefore not formally a Kitten state. There are two 
features of this state that are of particular interest. First, 
looking at the second term, it is evident that the photon 
pairs created by this operator are, in general, pairwise 
spectrally entangled. Second, in general the photon cre- 
ated by the first term has a different spectral form to the 
other photons. 

For the preparation of " real" Schrodinger Kitten states 
we are interessted in the conditions under which the en- 
tanglement between the different components vanishes. 
From the engineering of spontaneous parametric down- 
converion we know that spectral entanglement of bi- 
photon states can be eliminated by a careful design of 
the generation process Q. In this case the initial state 
from the down-conversion process is separable into iden- 
tical components (i.e. formally a superposition of Fock 
states). Thus Ay becomes separable, i.e. Xij = A^A'-, and 
it follows that the photon-subtracted state is of the form 



OO 



A" 



E^ 



2n-l 



|0>, 



(62) 



where the leading factor of Ao reflects the non- 
dctcrminism of the conditioning process. Clearly in this 
case the state is formally a superposition of Fock states, 
as required. Note, that in this second analysis we chose 
our mode decomposition according to the internal struc- 
ture of the parametric downconversion states. Then the 
mode of the photon which we use for post-selection gets 
decoupled from the signal state and narrow band filtering 
becomes obsolete. This provides us with a clear strategy 
to realize Schrodinger kitten states in practise such that 
the spatio-spectral structure of the states can be ignored 
and an effective single mode treatment becomes valid. 

However, in general the down-conversion process does 
not produce separable photon pairs, but highly spectrally 
entangled pairs as a result of energy conservation con- 
ditions. Thus, this technique will not formally produce 
Kitten states, but rather a more general class of states 
exhibiting the same photon number distribution as Kit- 
ten states. We might be tempted to utilize narrowband 
filitering for purifiying the state. Though, recall from the 
previous sub-sections that spectral filtering is equivalent 



to a decomposition into one signal and one unobserved 
eigenmode with a subsequent tracing over the unobserved 
mode. When a photon pair is spectrally entangled it is 
not possible for both photons to simultaneously have per- 
fect mode overlap with the defined signal mode - an en- 
tangled state cannot have unit overlap with a separable 
one. This will degrade the purity of the states with corre- 
lations between A\ , Aj . When performing homodync de- 
tection, imperfect mode overlap with the coherent probe 
beam also leads to a tracing out of the orthogonal compo- 
nents of the measured mode. This effectively introduces 
mixing with the vacuum state. Thus, in general we ex- 
pect this procedure to result in observed mixing with 
the vacuum state, even in the presence of perfect mode- 
matching, i.e. equivalent spectrum for sigal and local os- 
cillator mode. Note that in the experiment of Ref. [It} 
significant mixing of this type was observed, ranging from 
29%— 36%. Based on our results we expect that this arises 
from a combination of mode-mismatch and spectral en- 
tanglement effects. It is important to recognize that the 
later effect is, in general, inherent to this state prepa- 
ration technique. This type of mixing cannot be com- 
pensated for through filtering, post-selection or improved 
alignment and should be distinguished from experimental 
imperfections. 

In conclusion the fideltiy for the conditional prepara- 
tion of Schrodinger kitten states is already constrained 
by the generation of the squeezed states. Spectral cor- 
relations between differnt photon numbers degrade the 
achievable purity such that the best avenue for improv- 
ing state fidelity is to engineer superpostions of true Fock 
states in the first instance. 



E. Conditional preparation of Fock states via 
non-degenerate parametric down-conversion 

Finally we consider the example of conditional 
state preparation via non-degenerate parametric down- 
conversion. This technique is widely used in the prepa- 
ration of single photon states for quantum optics appli- 
cations, including quantum information processing ones. 
We consider the more general case where we postselect 
upon detecting n photons in one arm of a non-degenerate 
down-converter and examine the form of the conditional 
state in the other arm. In the non-degenerate case the 
Hamiltonian differs only slightly from Eq. [57J an d takes 
the form, 

H(t) = A / / ijiiuii, uj2)a < (lj\)w{(jJ2) doji &UJ2 + H.c. 



Following a similar derivation as per the previous exam- 
ple, the prepared state is of the form 



E 

71 = 



E A <-- 1 >\ 



|0>. 



(64) 
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Now suppose we post-select upon detecting m photons in 
mode A. In general the detection process may involve fil- 
tering. We will use the usual trick of choosing our eigen- 
mode basis to contain the detector mode as one of its 
elements, £o- Thus, the detector effectively applies the 
projector (|£o)(£o|) m - Upon applying this projector, the 
conditionally prepared state is 



IV. 



A" 



cond/ 



|0>, 



(65) 



where we ignore terms n > m. This approximation is 
valid as long as A <C 1, which is typically the case if the 
downconversion source is operated in the low gain regime. 
We can relabel the term in the bracket to obtain 



^rM |0) ' 

to! V / 



(66) 



where — ^ . XojBj is the mode creation operator 
characterized by the marginal distribution of the joint 
SDF over mode B. Importantly, in contrast to the previ- 
ous example the conditionally prepared state in Eq. 1661 is 
always formally an m-photon Fock state, independent of 
A making this approach to state preparation very suitable 
for quantum optics applications. Note that this feature 
arises from the particular form of Eq. 1641 specifically the 
pairwise spectral entanglement, whereas for a completely 
arbitrary multi-photon state this property will not apply. 
This highlights the variability in the behavior of systems 
under filtering and measurement, and the need to ana- 
lyze them on an individual basis. In some cases, as per 
the previous examples, filtering has a major impact on 
the form of the final state of the system, while in others, 
as here, it does not. 

In this analysis we have assumed the detector projects 
onto a particular, well-defined mode |£q)- As mentioned 
in the introduction realistic detectors with single photon 
sensitivity, namely avalanche photo diodes, lack the time 
resolution to distinguish the temporal modes. This leads 
to the annihilation of off-diagonal frequency modes Q- 
For this reason narrowband filtering is an effective tool 
in increasing the purity of conditionally prepared states, 
since it effectively forces the detector to respond to only 
a single well-defined mode [27| • Our result of this subsec- 
tions clarify that the the strategy to improve the purity 
by narrow band filtering, which is known for single pho- 
ton state, can be extended to Fock states of higher or- 
der number states if no superpostions of different photon 
numbers are considered. 



V. CONCLUSION 

We have discussed the spectral structure and decompo- 
sition of multi-photon states. We considered the general 
properties of the spectral distribution function and in- 
troduced an approach for decomposing it into discretized 



orthogonal 'eigenmodes'. This leads naturally to an occu- 
pation number representation for general optical states. 
The advantage of this representation is that it simplifies 
many types of calculations. There are two main reasons 
for this. First, we have great freedom in choosing the basis 
for the decomposition. So, for example, when consider- 
ing the interaction between two states, we can decompose 
one state into a basis consisting of components overlap- 
ping with and orthogonal to the other state. This allows 
for very straightforward treatment. Second, because the 
discrete eigenmodes are orthogonal, many types of evolu- 
tion (e.g. all linear optical interactions) can be calculated 
by acting the relevant evolution independently on each 
of the eigenmodes. This allows us to adopt an effective 
single-mode representation where we choose one under- 
lying spatio-spectral reference mode and define a photon 
number respresentation corresponding to this mode. We 
have shown that homodyne detection can be understood 
in this picture. Here mode mismatch is to be interpreted 
not only as an experimental imperfection but as a limita- 
tion inherent to a system comprising modes with corre- 
lated or different spatio-spectral characteristics. We have 
applied these techniques to several other examples, which 
provides further insight into the impact of the internal 
structure of multimode states on quantum state manip- 
ulation. For the Hong-Ou-Mandel interference our meth- 
ods highlight very clearly that the symmetry of states 
defines the visibility of the interference. Our analysis of 
spectral filtering as well as the Schrodinger Kitten state 
preparation proves that spectral filtering in general does 
not purify quantum states. The detection of vacuum mix- 
ing for conditionally prepared Schrodinger Kitten states, 
which is commonly attributed solely to experimental im- 
perfections, also arises from the spatio-spectral correla- 
tion of the internal structure of states. To eliminate these 
effects the best strategy is to engineer pure squeezed 
states which do not exhibit spatio-spectral correlations. 
In contrast to this example the generation of Fock states 
of higher photon number can be accomplished by narrow 
band filtering of a parametric downconversion state. 



In conclusion, our examples illustrate that different 
optical systems have to be analyzed on an individual 
basis. In many situations an adapted mode decompo- 
sition simplifies calculations and provides the means to 
identify the impact of spatio-spectral structure. We be- 
lieve the presented techniques will prove useful in the 
theoretical analysis of quantum optical systems, particu- 
larly in situations where 'non-ideal' effects such as mode- 
mismatch, photon distinguishability and imperfect detec- 
tors are considered. 
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APPENDIX A: PROOF THAT M = n\ IF AND 
ONLY IF ip(.) EXHIBITS FULL PERMUTATION 
SYMMETRY 



APPENDIX C: FOUR-PHOTON 
NORMALIZATION FACTOR IN FOUR-PHOTON 
INTERFERENCE 



By definition, 



i «n)' E ^(i > [w l ,...,Wn])d(3. (Al) 

-Pes„ 



We will adopt the shorthand ip = ip(u>i, . . . , w n ) and 
?/>p = ^(F[wi,...,w„]). We now employ the Cauchy- 
Schwarz inequality, 



K^p>| < IMI-lhMI, 



(A2) 



where H^H = (ip,ip)- In our case we have \\ip\\ = \\ipp\\ — 
1. Thus, 



\{1>,1>p)\ < 1, 



(A3) 



with equality for ip — ipp. The expansion for Af con- 
tains n! terms, each of the form (ip, ipp) which have mag- 
nitude less than or equal to unity. Furthermore, since 
N = (ip\ip),N e R. In general (ip,ip P ) G C, but if N = n\ 
then necessarily (ip, ip P ) = 1 V P, i.e ip = ipp V P. 



APPENDIX B: PROOF THAT A BASIS 
CONTAINING AN ARBITRARY INTEGRABLE 
FUNCTION AS ONE OF ITS ELEMENTS 
ALWAYS EXISTS 

Let Tt be a separable Hilbert space with countable in- 
finite basis X = {x n }. Let / be a normalized function 
TL — > TL. We wish to show that the set X' = X\xk U / 
is a linearly independent set which can be used to create 
the orthonormal basis Y = {y n } with yo = f. 

Since f £ TC and X forms a basis for 7i there must exist 
at least one Xk such that (xk, /) =/= 0. Then X\xk remains 
linearly independent but is no longer complete. Adding / 
to X\xk restores completeness, but X' = X\xkLif is not 
an orthonormal basis because there exist i ^ k such that 
(xi, /) 7^ 0. Because X is countable, the number of expan- 
sion coefficients for / is also countable, so the expansion 
of any g € TL is possible with a finite number of expan- 
sion coefficients using X' . The integral form of Holder's 
inequality (see, for example, Ref. [38|) guarantees the ex- 
istence of the scalar product defined over X' . Since X' 
is linearly independent, an orthonormal basis Y can be 
constructed using the Gram-Schmidt orthogonalization 
procedure. Note that this argument relies on the separa- 
bility of the underlying Hilbert space. This is fulfilled for 
most systems of interest, but there are counterexamples 
(the standard one being an infinite spin chain) which do 
not fulfill this requirement. In such cases, convergence 
issues will need to be investigated separately. 



We wish to calculate the normalization factor A/4 for a 
state of the form 

|V>4>= 4,4,1^4, |0>. (CI) 

We have 

AA 4 = (^4) = (OlA^A^A^A^A^A^AlJlJO). 

(C2) 

Summing over all combinations in which terms from the 
left hand side of the expansion can act on terms from the 
right hand side we obtain 



A/4 = 47n 2 722 2 + 4711722712721 + 47i2 2 72i 2 , 
where 

7tf = (014^10). 
Note that 711 = 722 = 1- Thus, 

A/4 = 4(1+7 + 7 2 ), 

where 



7 = 712721 



J tfi(uj)*ip2{uj) 



did 



(C3) 
(C4) 
(C5) 

(C6) 



APPENDIX D: DERIVATION OF THE PDC 
HAMILTONIAN 

Following [3^1 , the time-dependent Hamiltonian for the 
PDC process is given by 

H{t)oc^^ij(k s ,h,t)ai(k a )a}(h)+h.c. (Dl) 

where h.c. denotes hermitian conjugation and ip(k s , fc;, t) 
characterizes the joint spatio/temporal structure of the 
signal and idler photons. Collecting all creation operators 
in H('(t) and all annihilation operators in iP + )(i), we 
can write the Hamiltonian as 



ff(t) = #<->(*) + £<+>(t),. 



(D2) 



Time propagation of the output state is governed by 
the Schrodinger equation. We obtain a solution to the 
differential equation using the standard integral ansatz 
which can be found in any standard textbook on quantum 
mechanics (e.g., [i(|). This leads to state of the form 
IV>> = E,>n> = E„£n|0> where 



U'o) =|0) 



VPrx 



dri / dr 2 ■ 

—00 J — 00 



dr„ 



(D3) 



(D4) 
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Consider the product of Hamiltonians for different 
times Tj,: 



II V(Tk) = l[(H ( -HTk) + # + >fa)) (D5) 

(D6) 



fc=i 



fc=i 



= £ n^ (s[fc % fe ) 

se(±)* fc=i 



Assuming that the number of generated photons is neg- 
ligible compared to the number of possible field modes, 
annihilation operators contained in H ( + ) will predomi- 
nantly hit a vacuum mode. Thus nLi H^\T k )\Q) = 
if there is at least one k such that S[k] — "+". Only 
contributions with S[k] = "-" for all k will effectively 
contribute. This simplifies Eqn. ID4l as follows. 



dn / dr 2 • • • 

— oo J — oo •/ — X 

xffH( n )...ff(-)( Tfl ) 
M ' ' ' drff(~>(r) 



dr„ 



(D7) 
(D8) 



n=0 



dwidw2V'(wij ^2)0 (1^1)0^(^2, 



(D9) 



This is the form of the Hamiltonian conventionally 
employed for the analysis of parametric downconversion 
which we also employ in this paper. 
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